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High-throughput ab initio analysis of the 
Bi-ln, Bi-Mg, Bi-Sb, In-Mg, In-Sb, and Mg-Sb systems 
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Prediction and characterization of crystal structures of alloys are a key problem in materials re- 
search. Using high-throughput ab initio calculations we explore the low-temperature phase diagrams 
for the following systems: Bi-ln, Bi-Mg, Bi-Sb, In-Mg, In-Sb, and Mg-Sb. For the experimentally 
observed phases in these systems we provide information about their stability at low temperatures. 
Keywords: Binary Alloys, Ab initio, Intermetallics, Transition Metals, Structure Prediction, Phase 
Stability, Magnesium, Indium, Bismuth, Antimony. 
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I. INTRODUCTION 

Magnesium alloys present an excellent combination of 
light weight and specific strength compared to steel and 
aluminum alloys and would offer an opportunity for radi- 
cal improvement in automotive vehicle design and perfor- 
mance [1,2]. However, for more than fifty years magne- 
sium alloy research has been confined almost exclusively 
to casting alloys and there has been no equivalent effort in 
developing wrought magnesium alloys of sufficiently high 
strength [3] . Both wrought aluminum alloys and wrought 
magnesium alloys achieve their increased strength lev- 
els through the process of precipitation-hardening, also 
known as age-hardening. In the case of aluminum al- 
loys, the metallurgy of precipitation-hardening has been 
thoroughly explored and has resulted in the five and six 
element aluminum alloys which form the basis for all 
aerospace structures [4] . The trial and error research has 
required an enormous effort over many years to produce 
the aluminum alloys now in use. This is so because the 
hardening process depends on the development of stable 
and metastable phases in the base alloy matrix, a prob- 
lem particularly difficult for experimental investigation 
due to long annealing times involved in the phase for- 
mation. If the development of greatly increased strength 
in precipitation-hardenable magnesium alloys followed a 
similar trial and error path it would require a similar level 
of effort [5]. 

Fortunately, the combination of ab initio density func- 
tional theory methods and data mining techniques pro- 
vides an opportunity to dramatically accelerate materials 
research by efficiently predicting new phases and accu- 
rately describing their ground states [6-10]. These the- 
oretical methods are particularly suitable for investiga- 
tion of low-temperature compounds and can thus play 
an important role in the development of precipitation- 
hardenable magnesium alloys. 

In this paper we systematically explore the low- 
temperature phase diagrams for the following six binary 
alloy systems: Bi-ln, Bi-Mg, Bi-Sb, In-Mg, In-Sb, and 
Mg-Sb. Both Mg-Bi and Mg-Sb systems have the poten- 



tial to be new age-hardenable magnesium-based alloys. 
Indium is of interest because it is one of the very few 
elements with high solubility in magnesium and yet has 
low electrochemical activity; whereby the relatively high 
electrochemical activity of magnesium might reasonably 
be expected to be reduced, albeit at an increase in overall 
alloy density. The In-Sb and In-Bi systems are included 
to complete the possibilities for binary alloys among these 
elements. 

The Binary Alloy Phase Diagrams [11] and the Paul- 
ing File [12] give a broad review of experimental data 
on these systems. In most cases the experimental results 
are complemented by thermodynamic modeling [11-13]. 
However, to the best of our knowledge, ab initio stud- 
ies regarding phase stabilities in these systems are scarce 
[14-16]. To check the completeness of the experimentally 
known phase diagrams in these systems and provide in- 
formation on the phase stabilities at low temperatures we 
have chosen a large library of most common prototypes 
in binary alloys and calculated ground state energies of 
these structures with ab initio methods. 

We describe the prototype library in the next section, 
give the details of the ab initio methods used in this 
study in section "High-throughput First Principles cal- 
culations" , and present the results for each of the six 
systems in section "Alloys" . 



II. BINARY SYSTEMS AND STRUCTURE 
PROTOTYPES 

We calculate six alloys, Bi-ln, Bi-Mg, Bi-Sb, In-Mg, 
In-Sb, and Mg-Sb in 186 crystal structure configurations. 
Many of these configurations have the same prototype, 
for example, AB3 and A3B so the number of distinct 
prototypes is 108. The various concentrations are listed 
in the following table. 
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Compounds 


Cone. 


number of 


composition 


of B 


prototypes 


A & B 


0% 


6 


A 5 B & AB 5 


16.66% 


3 


A 4 B & AB 4 


20% 


2 


A 3 B & AB 3 


25% 


27 


A 2 B 5 & A 5 B 2 


28.57% 


1 


A 2 B & AB 2 


33.33% 


34 


A5B3 & A3B5 


37.5% 


3 


A3B2 & A 2 B 3 


40% 


2 


A4B3 & A3B4 


42.85% 


1 


AB (& BA*) 


50% 


29 (+3) 



TABLE 1. Compositions, concentrations and number of 
prototypes in the library. The library has 186 structures, 
and 108 distinct prototypes (* at composition AB, 3 pro- 
totypes have different point groups in atomic positions A 
and B, therefore they represent distinct structure types). 

Of such prototypes, 67 are chosen from the most com- 
mon intermetallic binary structures in the Pauling File 
[12] and the CRYSTMET database [13], plus the common 
low temperature compounds (5) reported for the Mg, In, 
Sb and Bi systems [11,12]. Such prototypes can be de- 
scribed by their Strukturbericht designation and/or nat- 
ural prototypes [11,12]: Al, A2, A3, A4, A6, A7, A15, 
B h , Bl, B2, B3, B4, B81, B8 2 , B10, Bll, B19, B27, B32, 
B33 (B/), C c , C2, C6, Cll 6 , C14, C15, C15 6 , C16, C18, 
C22, C32, C33, C37, C38, C49, D0 a , D0 3 , D0 9 , D0n, 
DO19, D0 22 , D0 23 , D0 24 , Dl 3 , Dl„, D2 d , D5 2 , D7 3 , D8 8 , 
D8 S , Lie, Lli, Ll 2 , L6 , Caln 2 , Cr 3 B 5 , CuTe, CuZr 2 , 
GdSi 2 (1.4), Mg 2 In, MoPt 2 , NbAs (NbP), NbPd 3 , Ni 2 In, 
Ni 2 Si, n (with z=l/4), Pu 3 Al (C03V), Ti 3 Cu 4 , W 5 Si 3 , 
YCd 3 , ZrSi 2 , 7-Ir. The rest of the structures (36) are fee, 
bec or hep superstructures. Twelve of these superstruc- 
tures consist of stacking of pure A and B planes along 
some common direction [6,17]. 



III. HIGH-THROUGHPUT FIRST PRINCIPLES 
CALCULATIONS 



were performed using at least 3500/ (number of atoms 
in unit cell) k-points distributed on a Monkhorst-Pack 
mesh [21,22]. With these energy cutoffs and k-points 
meshes the absolute energy is converged to better than 
10 meV/atom. Energy differences between structures are 
expected to be converged to much smaller tolerances. All 
structures were fully relaxed. 

PAW-GGA calculations. When several structures 
are in close competition for the ground state, we also per- 
formed calculations in the Generalized Gradient Approxi- 
mation (GGA), with Projector Augmented- Wave (PAW) 
pscudopotentials, as implemented in VASP [20,23,24]. In 
general, we expect the PAW-GGA approach to be more 
accurate than the US-LDA. For the GGA correlation en- 
ergy, we used the Perdew-Wang parameterization (GGA- 
PW91) [25]. Compared to the US-LDA case, we use an 
increased energy cutoff of 1.75 times the larger of the 
suggested energy cutoffs for the elements in the system 
and a finer k-point mesh with at least ^4000/ (number 
of atoms in unit cell). 

Symmetries of the pure elements. Our calcula- 
tions reproduce the correct experimental crystal struc- 
tures of the pure elements at room temperature. Bi and 
Sb are most stable in the A7 structure (a- As prototype), 
while Mg and In have hexagonal closed packed (A3) and 
face-centered tetragonal (A6) structures, respectively. 

Calculation of the formation energies and the 
convex hull. The formation energy for each structure 
is determined with respect to the most stable structure 
of the pure elements. To determine the ground states 
of a system one needs to find, as a function of composi- 
tion, the ordered compounds that have an energy lower 
than any other structure or any linear combination of 
structures that gives the proper composition. This set of 
ground state structures forms a convex hull, as all other 
structures have an energy that falls above the set of tic 
lines that connects the energy of the ground states. In 
thermodynamical terms, the convex hull represents the 
Gibbs free energy of the alloy at zero temperature. 



The high-throughput ab initio method used for this 
project is fully described in references [6,7] A summary 
of the details of the calculations is given below. 

Ultra Soft Pseudopotentials LDA calculations 
(US-LDA). The energy calculations were performed us- 
ing Density Functional Theory in the Local Density Ap- 
proximation (LDA), with the Ceperley- Alder form for the 
correlation energy as parameterized by Perdew-Zungcr 
[18] and with ultra soft Vandcrbilt type pscudopotentials 
[19], as implemented in VASP [20]. Calculations are done 
at zero temperature and pressure, with spin polarization, 
and without zero-point motion. The energy cutoff in an 
alloy was set to 1.5 times the larger of the suggested en- 
ergy cutoffs of the pseudopotentials of the elements of 
the alloy (suggested energy cutoffs are derived by the 
method described in [20]). Brillouin zone integrations 
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IV. ALLOYS 

Bi-In (Bismuth - Indium) 

The phase diagram of the Bi-In system is known from 
experimental investigations and thermodynamic model- 
ing [11,12,26-32]. Three phases are expected to be stable 
at low temperatures: Biln-BlO, BiIn2-InNi2, and Bi3±n5- 
Cr 3 B 5 . However, the evidence for the low-temperature 
stability of Biln2 and Bi3iri5 is not conclusive [27]. With 
our US-LDA calculations, we have found that only three 
phases have negative formation energy. These phases 
Biln-BlO, BiIn2-InNi2, and a monoclinic structure with 
space group C2/m #12 at 50% concentration have for- 
mation energies equal to -17 meV/atom, -4.4 meV/atom, 
and -7.9 meV/atom, respectively. Therefore, phase dia- 
gram calculated in US-LDA has only one stable com- 
pound, Biln-BlO; the other two experimentally observed 
stable phases Biln2 and Bialns are above the convex hull 
by 6.3 and 26 meV/atom, respectively. 



Bi-In system 


Low Temperature Phases comparison chart 


Composition 
% Bi 


Experimental 
(Massalski [11]) 


A b initio 
result 


33.3 


BiIn 2 -InNi 2 


two phase region (us-lda) 
BiIn2-InNi2 is 

~6.3 meV/at. higher than 
the tie-line In <->BiIn 
Biln2 stable (paw-gga) 


37.5 


BisIns-CrsBs 


Bi3lns unstable 
us-lda and paw-gga 


50 


Biln-BlO 


Biln-BlO 
E /= -17 meV/at. (us-lda) 
Ef=-6meV/at. (paw-gga) 




40 60 80 100 

In Atomic Percent Bismuth Bi 

FIG. 1. Bi-In (Bismuth - Indium) ground state convex hull 

calculated in US-LDA. 

Due to the rather small energy differences and overall 
low formation energies we further investigate the phase 
diagram with the PAW-GGA potentials, as described in 
the method section. With PAW-GGA, we find that Biln- 
BlO has less negative formation energy, -6 meV / atom, 
and therefore BiIn2-InNi2 becomes stable by having its 
energy 3.1 meV/atom below the tie-line In <-^>BiIn. The 



monoclinic structure with space group C2/m #12 at 
50% concentration is now unstable having the forma- 
tion energy of 10 meV/atom. The Bi^In^ phase still 
has positive formation energy (22 meV/atom) and re- 
mains significantly higher the tie-line InNiq <^InBi (by 
29 meV/atom). In summary, our results suggest that at 
low temperatures Biln 2 may be stable, while Bi 3 In 5 is 
likely to be unstable. 



0.01 
0.005 
c 

£ A6 

% -0.005 
> 

CD 

-0.01 
-0.015 
-0.02 









C2/m 




B10 


lnNl s 





20 40 60 80 100 

In Atomic Percent Bismuth Bi 

FIG. 2. Bi-In (Bismuth - Indium) ground state convex hull 

calculated in PAW-GGA. 



Bi-Mg (Bismuth - Magnesium) 

Only one stable compound has been observed for 
the Mg-Sb system at low temperatures: ctBi 3 Mg 2 - 
La3C>2 [11,12,33-35]. Our calculations confirm the low- 
temperature stability of this phase with a formation en- 
ergy of -252 meV/atom. We have not found any other 
stable compounds for this system, therefore the experi- 
mental low temperature part of the diagram is complete. 



Bi-Mg system 


Low Temperature Phases comparison chart 


Composition 
% Bi 


Experimental 
(Massalski [11]) 


Ab initio 
result 


40 


D5 2 


Bi 3 Mg 2 -D5 2 (us-lda) 
E /= -252 meV/at. 
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FIG. 3. Bi-Mg (Bismuth - Magnesium) ground state con- 
vex hull. 
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Bi-Sb (Bismuth - Antimony) 

No intermediate stable compounds have been reported 
for the Bi-Sb system and it is considered to be a non- 
compound forming system [11,12]. Our calculations con- 
firm the absence of any stable phases. The structure with 
lowest formation energy in the whole range of concentra- 
tions is found to be BiSb-Bl with E/=21 meV/atom. 

In-Mg (Indium - Magnesium) 

Five ordered compounds have been reported for 
the system Mg-In at low temperatures [11,12,36-41]: 
M gl . 2 In 2 . 8 -Ll 2 (7'), Mgln-Llo (/?"), Mg 2 In (prototype, 
ft, Mg 2 Tl in [11]), Mg 5 In 2 -D8 g (ft), and Mg 3 In (ft) 
with Pearson symbol hR48, and space group R3m (# 
166) [11,12] (note that ft has Pearson symbol hR16 in 
Massalski [11]). At high temperatures phase ft trans- 
forms into the Ll 2 structure, which is similar to the high 
temperature field ft. Phases ft and ft share a eutectoid 
reaction at 202°C (ft «-> ft + ft at ~ 27.5% of In), and 
a perictectoid reaction at 337° C (ft <-> (3 + ft at ~ 24% 
of In). Therefore we expect the morphology of ft to be 
similar to Ll 2 (with partial occupation) and to find the 
structure Mg3in-Ll 2 under the tie-line I11-A6 <->MgIn- 
LIq at zero temperature. 



In-Mg system 



Low Temperature Phases comparison chart 



Composition 
% In 


Experimental 
(Massalski [11]) 


Ab initio 
result 


~26 to 38.5 


Mg 3 In (ft) 
hR48 R3m [12] (low T), 
~ Ll 2 (high T) 


Ll 2 


28.6 


Mg 5 In 2 -D8 9 (/3 3 ) 


two phase region 
D8 9 is 4.9 meV/at. 
above ft *-> ft (us-lda) 

D8 g is 5.4 meV/at. 
above ft <-» ft (paw-gga) 


~34 


Mg 2 In (ft) 


Mg2ln stable 
~3.5 meV/at. below 
LI2 <->Llrj (us-lda) 
~1 meV/at. below 
LI2 *->L^o (v aw ~99 a ) 


~39 to 59 


Mgln-Llo (13") 


Ll 


-69.5 to 75.5 


Mgi.2ln2.8-Ll2 (7') 


Mgln 3 -Ll 2 



The off-stoichiometry Mgi. 2 In 2 .8-Ll 2 phase is not sub- 
ject of our investigation, since it requires simulations 
of disordered systems. With our calculations, we find 
the ordered phase Mgln3-Ll 2 to be stable. Two com- 
peting Long Period Superstructures of Mgln 3 -Ll 2 , D0 2 3 
and D0 2 4, are less favorable than the Ll 2 phase by ^2.5 
and ~5.7 meV/atom, respectively. This suggests that 
long-range interactions are weak in this system at com- 
position Mgln3. We confirm the stability of Mgln-Llo 
(ft'), Mg 2 In-ft, and Mgln 3 -Ll 2 phases but find that 
Mg5ln 2 -D8 s lies 4.9 meV/atom above the tie-line Mg3ln- 
ft <-^Mg 2 In-ft . Besides, Mg 2 In-ft is stable but with a 
relative stability energy of 3.5 meV/atom below the tie- 
line Mg3ln ^->Mgin, which is within the error of present 



calculations. We find a metastable compound, Mgln 2 -C c , 
to be ~2.7 meV/atom above the tie-line Mgln <->MgIn3. 
This compound might become stable at higher temper- 
atures and pressures. We recalculate the phase diagram 
with PAW-GGA potentials and find essentially no dif- 
ference compared to the US-LDA results: M<7 2 /n-ft is 1 
meV/atom below the Mg 3 In <^>MgIn tie-line, Mg 5 In2~D8 g 
is above the Mg^In-j3\ <->M<? 2 /n-ft by 5.4 meV/atom and 
Mgln2-C c remains metastable by 4.6 meV/atom. In sum- 
mary, our calculations confirm the low temperature part 
of the Mg-In diagram known from experiment. The re- 
sults show that Mg5ln2-D8 g may not be stable at low tem- 
peratures. In addition, we identify a metastable phase, 
Mgln2-C c , which might be stable at higher pressures and 
temperatures. 



0.02 




20 40 60 80 100 

Mg Atomic Percent Indium In 

FIG. 4. In-Mg (Indium - Magnesium) ground state convex 

hull. 
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In-Sb (Indium - Antimony) 

The phase diagram of the system In-Sb is based on a sin- 
gle compound aInSb-B3 [11,12,42-47]. It is well studied 
at high pressures and temperatures where other modi- 
fications of InSb (/3InSb, 7lnSb, and <5InSb) have been 
observed [28,45]. The high-pressure transitions have also 
been thoroughly studied with ab initio methods [15,16]. 
At zero pressure our calculations confirm the experimen- 
tal stability of the low temperature compound odnSb-B3. 
The InSb-BlO phase, stable for the similar binary system 
In-Bi at 50% concentration, here is unstable with an en- 
ergy ~105 meV/atom higher than aInSb-B3. We do not 
observe any other stable compounds for the In-Sb sys- 
tem and conclude that the low temperature experimental 
characterization of the system is complete. 



In-Sb system 


Low Temperature Phases comparison chart 


Composition 
% Sb 


Experimental 
(Massalski [11]) 


A b initio 
result 


50 


B3 


aInSb-B3 
B10~105 meV/at. 
higher than B3 (us-lda) 
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FIG. 5. In-Sb (Indium - Antimony) ground state convex 
hull. 

Mg-Sb (Magnesium - Antimony) 

Similarly to the Bi-Mg system the Mg-Sb system 
has only one stable compound at low temperature: 
aMg 3 Sb 2 -B5 2 [12,11,48,35,49-51]. Interestingly, there 
is no solid solubility for Mg in Sb or Sb in Mg, al- 
though Mg and Sb have similar dimensions (atomic ra- 
dius ratio Sb/Mg=l. 09) [11]. At high temperatures phase 
aMg3Sb2-B52 undergoes a polymorphic change, as re- 
ported in [35]. To our knowledge, there is no information 
about the stability of aMg 3 Sb 2 -B5 2 below 450°C and 
it is unclear how far this stable phase extends into the 
low temperature region. Our calculations confirm that 
this compound remains stable at low temperatures hav- 
ing the formation energy of -404 meV/atom. No other 



stable compounds have been found for this system, there- 
fore the experimentally known low-temperature part of 
the phase diagram is complete. 



Mg-In system 


Low Temperature Phases comparison chart 


Composition 
% Sb 


Experimental 
(Massalski [11]) 


Ab initio 
result 


40 


D5 2 


Mg 3 Sb 2 -D52 
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FIG. 6. Mg-Sb (Magnesium - Antimony) ground state con- 
vex hull. 



V. CONCLUSIONS 

The results of our systematic ab initio study of the 
six binary systems Bi-In, Bi-Mg, Bi-Sb, In-Mg, In-Sb, 
and Mg-Sb suggest that the experimentally known low- 
temperature phase diagrams for these alloys are com- 
plete. Using the calculated ground state energies we 
complement experimental information on the stability of 
several phases in the limit of low temperature and zero 
pressure: phase aMg 3 Sb 2 -B5 2 is stable, phase Mg 5 in 2 - 
D8 S may not be stable and phase Bialns is likely to be 
unstable. We find that compound Mgln 2 -C c , not ob- 
served experimentally, is metastable but very close to 
stability and could therefore exist at different tempera- 
tures and pressures. Further theoretical investigation of 
these structures and search for other (meta) stable phases 
at finite temperatures and pressures should include ther- 
modynamic effects. 
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